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Abstract 

The best quadratic approximation to the retarded polaron action due to Adamowski et at and 
Saitoh is investigated numerically for a wide range of coupling constants. The non-linear variational 
equations are solved iteratively with an efficient method in order to obtain the ground-state energy 
and the effective mass of the polaron. The virial theorem and expansions for small and large 
couplings are used to check the high accuracy of the numerical results. Only small improvements 
over Feynman's (non-optimal) results are observed. For a moving polaron it is shown that the 
most general quadratic trial action may contain anisotropic terms which, however, do not lead to 
improvements for the ground-state energy and effective mass. 
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1. The polaron problem is a non-relativistic field theory for an electron moving in a crystal and has 
received a lot of attention in the past decades (for reviews see [1], [2]). Among the many theoretical 
treatments Feynman's approach [3] is still outstanding: he first integrated out the phonons to obtain 
an effective action for the electron which for large Euclidean times reads 

s eS = / dt-i 2 --^ dtdt' — = — — . (i) 

Jo 2 2^/2 Jo | x(t) - x(t') | V ' 

He then performed a variational approximation with a quadratic retarded trial action 

S t = \ f dt x 2 + \ f dt dt' fit- [ x(t) - x(f ) ] 2 , (2) 

JO * JO 

choosing /f(ct) = C exp(— w), where C and it; are two variational parameters (the strength parameter 
C is usually written as w{v 2 — w 2 )/4 ). This yields one of the best analytical approximations for the 
ground-state energy of the polaron for all values of the dimensionless coupling constant a. The 
best possible, rotationally invariant, quadratic trial action is obtained by replacing the exponential 
retardation by an arbitrary function f(a) and was proposed by Adamowski et al. [4] and Saitoh 
[5]. Surprisingly for both small and strong coupling this best isotropic quadratic approximation only 
yields very small improvements for the ground-state energy Eq 

E 2=3 -a -a a 2 + 0(a 3 ) (3) 
-da 2 + 0(l). (4) 

One finds at small coupling ajr = 0.0123457, a- lso = 0.0125978 whereas the exact value is a = 0.0159196 
and for strong coupling ap = d\ so = 0.106103 compared to a = 0.108513 . Probably discouraged 
by these analytical results the best isotropic quadratic approximation has never been investigated 
systematically for a whole range of couplings, in particular for intermediate coupling. 

It is the purpose of the present Letter to do this for the ground-state energy and the effective mass 
of the polaron and to point out that an easy, efficient method exists to solve the non-linear variational 
equations. Although numerical results have been reported in the literature [4, 6], it will turn out that 
they are unreliable and considerably overestimate the improvement on Feynman's approach. We also 
investigate the question whether the inclusion of anisotropic terms in the quadratic trial action leads 
to further improvements. This work is an outgrow of recent attempts to generalize Feynman's polaron 
approach to four-dimensional field theories [7-11] in the context of the worldline formalism [12]. The 
nomenclature is the one used in Ref. [7]. 



2. We start with the expression for the ground-state energy using the quadratic trial action (2): 
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is the "pseudo-time" corresponding to the "profile function" A{E) which is linked to to the retardation 
function through 

A(E) = l + ^jf d*f(a)sm 2 (— J . (7) 

With Feynman's choice one obtains the standard expression Ap(E) = (v 2 + E 2 )/(w 2 + E 2 ) for the 
profile function and the pseudotime can also be given analytically. However, one can do better by not 
imposing a special form for the retardation function [4, 5]. Indeed, by varying Eq. (5) with respect 
to f(a) or, equivalently, with respect to A(E), one finds that the best variational profile function is 
determined by 

A iso (E) = l + / da^— K -=, T J - L . 8 

3v^ A) tf so (a) E 2 

By comparing with Eq. (7) one sees that the retardation function has the form 1 

/ iso (a) = -fl- . (9) 

In particular, since for small a the pseudotime [i 2 (cr) behaves like a, we see that the retardation 
function /i so (c) has a u _3 / 2 -singularity at small relative times in marked contrast to Feynman's 
Ansatz. Equivalently, the profile function at large E does not approach unity like 1/E 2 as Feynman's 
parametrization suggests but slower, like Ai SO (E) — > 1 + 2\/2a/(3E 3 ^ 2 ) . The wrong small-time 
behaviour of Feynman's Ansatz is responsible for the awkward behaviour of his variational parameters 
for a — > 0: both v and w tend to 3 in this limit whereas one would have expected w — ► 1 by comparison 
with the exact effective action (1). 

For a quantitative description at arbitrary a one has to solve the coupled non-linear variational 
equations (8) and (6) numerically. This can be done as follows [8]: one maps the infinite intervals 
to finite ones by substituting E = tan 2 , a = tan 2 ip and solves the variational equations on a 
grid of Gaussian points. In this way the required integrals both in the variational equations and in 
the evaluation of the ground-state energy can be evaluated directly by Gauss-Legendre quadrature. 
The above choice of mapping is made to eliminate the integrable but numerically slowly convergent 
square-root behaviour near a — ► in the various cr-integrals. The variational equations are then solved 
iteratively by starting either from the perturbative values A(E) = 1, /U 2 (cr) = a or from Feynman's 
parametrization using the known values of the variational parameters [13]. For larger coupling the 
latter method gives faster initial convergence as seen in Fig. 1. The convergence rate was monitored by 
calculating the maximal relative deviation between two iterations, both for A(E) and for n 2 (cr) for all 
Gaussian points, i.e. for all discrete values 0j and ipi. In the numerical calculations reported in Table 1 
the allowed maximal relative deviation was set to 10 -7 and reached after 7-70 iterations, depending 
on the coupling constant (see Fig. 1). Furthermore, the [0, 7r/2] -range for the variables 9,ip was 
subdivided into n e intervals and Gaussian integration with n g points was applied in the subintervals. 
Typically values of (n e ,n g ) = (6,72) were needed to obtain a ground-state energy accurate to six 
digits. It turned out that solving the non-linear variational equations by this method is numerically 
not more demanding than minimizing the energy functional with respect to the Feynman parameters 
v and w. 



1 It should be kept in mind that for the ground-state energy we only consider the (3 — » oo limit which leads to the 
already simplified expression given in Eq. (1). 
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Figure 1: Convergence rate for solving the variational equations (8, 6) iteratively at two different 
values of the coupling constant. The open symbols denote the case when the starting values are the 
free ones, the full symbols when Feynman's parametrization is used for initialization. 

To check the numerical stability of the method the ground-state energy was also calculated in a 
different way: the variational equations for A\ so (E) can be used to re-express the kinetic term in 
terms of the potential V (see the Appendix of Ref. [8] where this is detailed for the relativistic case): 



This can be combined with the potential term and after elimination of the derivative of the pseudotime 
one obtains the expression 



This may be called "virial energy" since it also arises from evaluating the virial theorem for two-time 
actions [14]. In the numerical calculations an agreement in six digits between the two expressions for 
the energy was demanded and obtained. 

Table 1 collects the results for coupling constants from a = 1 to a = 15. Note that convergence 
is obtained in all cases whereas in the relativistic case for coupling constants larger than a critical 
coupling the relative maximal deviations started to grow again with increasing number of iterations 
[8]. This was a signal for the instability of the scalar Wick-Cutkosky model whereas the polaron has 
a well-known strong coupling limit [15]. Also included is the effective mass of the polaron, defined by 
the expansion E(p) = Eq + p 2 /(2m*) + • • • . It is given by [6] 




(10) 




(11) 




(12) 



and coincides with the value of the profile function Ai so (E) at E = (see Eq. (8)). 

Being a variational calculation the values for the ground-state energy should lie below Feynman's 
result which is non-optimal. Table 1 shows that for small coupling the relative change of the energy 
follows the behaviour (of — ai so )a = —2.52 • 10~ 4 a which is expected from the small-coupling 
expansion (3). But also for intermediate coupling the best isotropic quadratic approximation gives 
only a slight improvement, maximally 0.15 % , for the ground-state energy. 
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-1.013296 


0.26 




1.19615 


0.53 
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-2.056467 


0.54 




1.47515 


2.21 




3 


-3.135951 


0.84 




1.89862 


5.12 




4 


-4.261309 


1.13 




2.60234 


8.93 




5 


-5.447781 


1.40 




3.93259 


12.1 




6 


-6.72130 


1.55 




6.9081 


10.2 




7 


-8.12440 


1.44 




14.416 


1.52 




8 


-9.70625 


1.12 


-0.70 


31.479 


-2.86 


-0.09 


9 


-11.49505 


0.81 


-0.75 


62.597 


-2.46 


-0.15 


10 


-13.4982 


0.58 


-0.78 


111.65 


-1.48 


-0.17 


11 


-15.7163 


0.41 


-0.79 


182.96 


-0.90 


-0.16 


12 


-18.1489 


0.30 


-0.79 


281.47 


-0.54 


-0.15 


13 


-20.7954 


0.23 


-0.80 


412.65 


-0.32 


-0.13 


14 


-23.6554 


0.17 


-0.81 


582.47 


-0.20 


-0.11 


15 


-26.7285 


0.13 


-0.81 


797.4 


-0.12 


-0.10 



Table 1: The polaron ground state energy and effective mass in the best isotropic quadratic approx- 
imation as a function of the coupling constant a. The third and sixth column display the relative 
deviation from the results using Feynman's parametrization (see Ref. [13] and Table 2), while the 
entries in the fourth and last columns should approach a constant value in the large coupling limit 
(see text). 

In contrast, the numerical values for this approximation reported by Adamowski et al. claim to be 
below the Feynman energy by as much as 0.75 % at a = ll. 2 That this is incorrect and probably due 
to insufficient convergence and/or numerical instability is shown in Fig. 2 where the present results 
for this specific coupling constant are plotted as function of the number of Gaussian points. It is 
clearly seen that with sufficient subdivisions of the intervals both expressions (5) and (11) for the 
ground-state energy converge to the same value which is far away from the one given in Ref. [4]. 

2 The difference might not seem much, however it should be remembered that Feynman's energy is itself only a little 
more than 2 % above the exact energy at asymptotically large couplings; see Eq. (4) . 
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Figure 2: Ground-state energy of the polaron for a = 11 as a function of the number n e of subdivisions 
used in the 72-point Gauss-Legendre integrations. Open circles denote the energy evaluated from Eq. 
(5), full circles from the virial theorem (11). The inset shows the convergence of both numerical values 
in more detail. The value using Feynman's parametrization and the result reported by Adamowski et 
al. [4] are shown by dashed lines. 



An additional check on the accuracy of the numerical results is provided by the large coupling 
limits of Eq° and Eq, which may be calculated analytically from Eq. (5). A systematic expansion 
around a — > oo may be developed by making use of the fact that as the coupling constant increases 
the pseudotime remains constant over an increasing range of a (see Sec. 3.2 as well as the Appendix 
of Ref. [9]). One finds, after some effort, that not only the leading [0(a 2 )] term of Eq° and Eq is 
identical (see Eq. (4)), but also the first subleading term: 

= ~<? " 3 Q + log 2) + O (a- 2 ) . (13) 

The numerical values in the fourth column of Table 1 indeed show that the difference between Eq° 
and Eq is of order oT 2 . 

For the effective mass the corrections to Feynman's result are somewhat larger (for small a the 
relative change follows the perturbative result 2(ai so — of) a 2 = 5.04 • 1CT 4 a 2 , and they are closer to 
the ones given in Ref. [6]. For example, at a = 5 (their largest value) Gerlach et al. give mf so = 3.940 
whereas a similar increase in Gaussian points as shown in Fig. 2 stabilizes the correct value at 
m* so = 3.93259. At larger couplings the convergence of the present method deteriorates somewhat 
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(see Fig. 1) but up to a = 15 four to five significant figures can be given also for the effective mass 
of Eq. (12). Surprisingly, the large coupling limits of both m* so and m F again coincide to the first 
subleading term, i.e., 

m^iso = « 4 - 3^ (1 + ^g 4) a 2 + O (a ) , (14) 

the numerical values quoted in the last column of Table 1 being consistent with this. In the course 
of these calculations it turned out that the values of the effective mass in Feynman's parametrization 
are very sensitive to the precise numbers for the parameters v, w and deviate somewhat from the ones 
given in Table II of Ref. [13]. We have re-calculated them by solving the variational equations for v, w 
(instead of minimizing the energy functional which is sufficient for the energy but not for the effective 
mass) and checked that they agree to a very high precision with the large-coupling expansion. For 
convenience and further reference the new values are given in Table 2. 
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3.10850 


2.86958 


-1.0130308 




1.19552 
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3.24535 


2.72564 


-2.0553560 




1.47189 
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3.42130 


2.56031 


-3.1333335 




1.88895 
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3.66464 


2.36792 


-4.2564809 




2.57931 
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4.03434 


2.14002 


-5.4401445 




3.88562 
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4.66687 


1.87363 


-6.710871 




6.83836 
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5.80989 


1.60365 


-8.112688 




14.3941 
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7.58682 


1.40329 


-9.695371 




31.5693 
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9.85025 


1.28230 


-11.485786 




62.7515 




10 


12.4749 


1.20918 


-13.490437 


-13.489166 


111.816 


111.150 


11 


15.4132 


1.16209 


-15.709808 


-15.709373 


183.125 


182.692 


12 


18.6483 


1.12988 


-18.143395 


-18.143230 


281.622 


281.327 


13 


22.1733 


1.10676 


-20.790681 


-20.790613 


412.782 


412.573 


14 


25.98515 


1.08952 


-23.651278 


-23.651248 


582.584 


582.432 


15 


30.08224 


1.076285 


-26.724904 


-26.724900 


797.498 


797.385 



Table 2: Precise values for Feynman parameters, ground state energy and effective mass. The values 
in columns 5 and 7 are from the strong coupling (SC) expansions E F (SC) = J2i=o e i % ■> m* F (SC) = 
J2i=o m i °^~ % f° r which we have calculated the following expansion coefficients: eo = —0.1061033, e 2 = 
-2.829442 , e 4 = -4.86387 ,e 6 = -34.1953 ,e 8 = 533.141 ,e 10 = 51525.2 ,e 12 = 6.61224 • 10 6 ,e 14 = 
9.53627 • 10 8 and m = 0.02001406 , m 2 = -1.012775 ,m 4 = 11.8558 ,m 6 = 43.0986 . Most of these 
coefficients, except e\ 2 and eu, have also been obtained in Ref. [16]. 
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3. For a moving polaron Eq. (2) is not the most general quadratic trial action. Indeed, there is 
now a preferred direction p = p/|p| which may be used for constructing a trial action which makes 
use of this directionality, i.e., 

Samso = S iso -i(X'-l)p- j o dt±(t) + ^J o dtdt'g(t-t')[p- [x(t)-x(t')]) 2 . (15) 

This the most general quadratic trial action which is scalar and invariant under time and space trans- 
lations as well as time-reversal. Time-translational invariance requires that the retardation functions 
in the quadratic terms are functions of the relative time and allows only a constant parameter (writ- 
ten as A' — 1) in the linear term. Space-translational invariance leads to a dependence on co-ordinate 
differences or derivatives 3 . Finally, time-reversal invariance requires the linear term to be purely 
imaginary. In Eq. (15) it has been written in such a way that it may be combined with the exponen- 
tial in the Fourier transform of the partition function projected on momentum p [14]. For p/Owe 
thus have a new variational parameter A' and two profile functions at our disposal 

A(E) — ► Aij(E) = A L (E)piPj + A T (E) (<% i,j = 1,2,3 . (16) 

Here the subscripts "L" and "T" denote longitudinal and transverse components with respect to the 
direction of the polaron momentum. Using the methods developed in Ref. [17] it is straightforward to 
calculate the various averages needed in the Feynman variational principle and to evaluate the limit 
[3 — > oo. One finds that the energy 

^aniso = l( 2 f2[A T ]+n[A L ])+Fp[AiL/4] + Y(2A-A 2 ) (17) 

is stationary 4 under variation of the parameter A = X'/Al(0) and the two profile functions Al^t(E) 
or, equivalently, the two pseudotimes pL 2 LT {a). The potential term is given by 

v = r dae -a r 1 ^ MgO 

P Jo ° 6 Jo X n\{?) + [v-t(o-) ~ x2 

Clearly, the symmetric solution At(E) = Al{E) = Ai so (E) and hence n\{a) = /u|(cr), emerges when 
p 2 = and hence, not surprisingly, E'Q iaso = Eq°. A little less obvious is that the effective mass will, 
in fact, also remain unchanged. To see this, we expand the exponential in the potential to order p 2 : 

\ 2 2 

V P [A T ,A L ,X] = V [A T ,A L ] + -JLvMtiAl] + O (p 4 ) . (19) 

To leading order in the momentum 

Vi[A r ,A L ] = Fi[A iso ] + o(p 2 ) = l-mf so + 0(p 2 ) (20) 

where Eq. (12) has been used. On the other hand, Vq (and f2) contain an additional term of 
order p 2 because of the implicit dependence of At,l on the momentum. By writing At,l(E) = 

3 Terms with derivatives may be converted into the form of Eq. (15) by suitable integrations by part. 

4 This is due to the use of a complex trial function and "momentum averaging" but as in Ref. [7] one can show that 
the energy for arbitrary momentum p is also minimal. The effective mass, however, is the derivative of the energy with 
respect to p 2 at p = and thereby has no minimum property anymore. 



(18) 
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A\ so (E) + p 2 AAt,l{E) + ... and expanding the functional, one can make make this remaining 
momentum dependence explicit. Hence one obtains for the energy 



iTianiso 
E P 



= n[A iso ] + V [A iso ] + P M dE 

Jo 



2 -^A T {E) + l -^A L {E) 



5A(E) 



n[A\ 



i=L,T 



- E + p 



5 ME) 
dE 



V [A T , A L }\ + ^ (2A - A 2 + AV^so]) + O (p 4 ) 

) A.— -Also 



2 r^ 2A ME) + AA L (E) 



5A{E) 



(n[A] + Vo[A 



A=A\. 



+ E_( 2 A-A 2 + A 2 y 1 [Aso]) + o( P 4 ) 



(21) 



0.000 



Wr* ★ ★ ★ ★ * X 5 

300 ° o o 



-0.005-^ 



Cl 



-0.010 - 



< 



-0.015 - 



-0.020 



J i I i I i L 



a 

•0.1 
■0.5 

o2. 

*5. 



o o 

second order: 
"■^expanded 



full \ 



-. — I — i — | — i — | — i — | — r- 

0.0 0.1 0.2 0.3 0.4 0. 

p 2 /2 



Figure 3: Energy difference between the variational approximation with the most general (anisotropic) 
and the isotropic trial action as function of the coupling constant a and the polaron momentum p. 
The numerical results at a = 5 have been multiplied by a factor 5 to make them better visible in the 
graph. The curves show the analytical perturbative 0(a 2 ) result: the dashed curve includes only the 
leading term at small momenta (i.e. 0(p 4 ), see Eq. (23) ) while the dotted line is the full second-order 
result without expansion in p 2 . 



The integral over E vanishes for any AAT t i(E) because of the variational equation for Ai so (E). On 
the other hand, the variational equation for A yields A -1 = 1 — ViL4i so ] so that, upon substitution of 
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Eq. (20), one obtains 



2m?. 



+ 



°(p 4 ) ■ 



(22) 



It is only at higher orders in p 2 that a difference between j^ 11180 and E^ can show up. In fact, by 
solving the variational equations perturbatively, one sees that there is a difference between these two 
energies at 0(p 4 ): 



T-iamso ttiiso 
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(23) 



The anisotropic trial action has lowered the energy, as it should, however note that numerically the 
improvement is very small, the dimensionless coefficient of the a 2 p 4 /4 term being only 0.004651. Fig. 3 
shows that for small couplings our numerical results obtained by solving the anisotropic variational 
equations indeed follow the analytical prediction (23) rather well. 



4. In summary, we have investigated the use of the most general quadratic trial actions, both isotropic 
and anisotropic, in Feynman's variational approach to the polaron problem and found out that they 
only lead to small numerical improvements compared to Feynman's original parametrization of the 
retardation function. This is despite the fact that the variational retardation functions have quite 
a different small-time behaviour compared to Feynman's Ansatz, but is consistent with results from 
the large- N expansion [18], where N is the number of space dimensions (similar results have been 
reported in Refs. [19]). Anisotropic terms in the trial action only become beneficial for a finite polaron 
momentum but not for the ground-state energy at p = and its derivative, the effective mass. 

The small improvement over Feynman's results, even with arbitrary retardation functions, indi- 
cates that a quadratic trial action leaves out part of the correct physics of the polaron. This is most 
evident when the second-order expression for the ground-state energy is compared with the exact re- 
sult. However, this does not imply that the optimal quadratic trial actions, both the isotropic and the 
anisotropic one, are useless. First, in four-dimensional relativistic field theory the correct ultraviolet 
behaviour of the trial action is of much greater importance [8] than for the polaron. Second, as shown in 
Ref. [9], the correct analytic structure of the profile function A(E) is essential for describing scattering 
processes, where an analytic continuation to normal time is required. For these purposes Feynman's 
Ansatz is totally inadequate whereas an improved parameterization fi(cr) oc exp(— w'a)/a 3 ^ 2 captures 
the essence of the variational retardation function (9) and allows an easy analytic continuation. Third, 
extensions to fermionic theories like Quantum Electrodynamics require precise but subtle relations 
between fermionic and bosonic retardation functions [11] which are only fulfilled by the variational 
solutions if the variational Ansatz is left sufficiently general. And, finally, for a massive relativistic 
particle on its gerlachmass-shell which corresponds to a moving 4-D polaron with p 2 = M 2 , the 4- 
momentum is large and therefore one may expect substantial improvements in describing the physics 
of the system (for example, its instability) if one uses the most general, anisotropic trial action. 



9 



Acknowledgement: One of the authors (R. R.) would like to thank M. Smondyrev for helpful 
discussions and remarks and G. Ganbold for pointing out additional references. 



References 

[1] B. Gerlach and H. L6wen, Rev. Mod. Phys. 63 (1991) 63. 

[2] T. K. Mitra, A. Chatterjee and S. Mukhopadhyay, Phys. Rep. 153 (1987) 91. 

[3] R. P. Feynman, Phys. Rev. 97 (1955) 660. 

[4] J. Adamowski, B. Gerlach and H. Leschke, in: J. P. Antoine and E. Tirapegui (Eds.), Functional 
Integration - Theory and Application, Plenum, New York, 1980, p. 291 . 

[5] M. Saitoh, J. Phys. Soc. Jpn. 49 (1980) 878. 

[6] B. Gerlach, H. Lowen and H. Schliffke, Phys. Rev. B36 (1987) 6320. 

[7] R. Rosenfelder and A. W. Schrciber, Phys. Rev. D53 (1996) 3337. 

[8] R. Rosenfelder and A. W. Schreiber, Phys. Rev. D53 (1996) 3354. 

[9] A. W. Schreiber, R. Rosenfelder and C. Alexandrou, Int. J. Mod. Phys. E5 (1996) 681. 

[10] A. W. Schreiber and R. Rosenfelder, Nucl. Phys. A601 (1996) 397; C. Alexandrou, R. Rosenfelder 
and A. W. Schreiber, Nucl. Phys. A628 (1998) 427; N. Fettes and R. Rosenfelder, Few-Body Syst. 
24 (1998) 1 . 

[11] C. Alexandrou, R. Rosenfelder and A. W. Schreiber, Phys. Rev. A59 (1999), 1762; Phys. Rev. 
D62 (2000) 085009 . 

[12] M. J. Strassler, Nucl. Phys. B385 (1992) 145; M. Reuter, M. G. Schmidt and Ch. Schubert, Ann. 
Phys. (N. Y. ) 259 (1997) 313 . 

[13] Y. Lu and R. Rosenfelder, Phys. Rev. B46 (1992) 5211 . 

[14] C. Alexandrou and R. Rosenfelder, Phys. Rep. 215 (1992) 1. 

[15] S. I. Pekar, " Untersuchungen zur Elektronentheorie der Kristalle", Akademie-Verlag, Berlin, 1954; 
S. Miyake, J. Phys. Soc. Jpn. 38 (1975) 181. 

[16] O. V. Selyugin and M. A. Smondyrev, Phys. Stat. Sol. (B) 155 (1989) 155. 

[17] R. Rosenfelder, C. Alexandrou and A. W. Schreiber, in: A. W. Schreiber, A. G. Williams and A. 
W. Thomas (Eds.), Proceedings of the Workshop on Methods in Non-Perturbative Field Theory, 
Adelaide (Australia), 2 -13 February 1998, World Scientific, Singapore, 1998, p. 163. 

10 



[18] M. A. Smondyrev, Physica A171 (1991) 191; M. A. Smondyrev, in: H. Grabert, A. Inomata, 
L. Schulman and U. Weiss (Eds.), Proceedings of the Fourth International Conference on Path 
Integrals from meV to MeV: Tutzing '92, World Scientific, Singapore, 1993, p. 190. 

[19] G. Ganbold and G. V. Efimov, Phys. Rev. B50 (1994) 3733; J. Phys. CIO (1998) 4845. 



11 



